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Abstract 

We describe the bound state and scattering properties of a quantum mechanical parti- 
cle in a scalar A-prong potential. Such a study is of special interest since these situations 
are intermediate between one and two dimensions. The energy levels for the special case 
of N identical prongs exhibit an alternating pattern of non-degeneracy and {N — 1) fold 
degeneracy. It is shown that the techniques of supersymmetric quantum mechanics can be 
used to generate new solutions. Solutions for prongs of arbitrary lengths are developed. 
Discussions of tunneling in A-well potentials and of scattering for piecewise constant po- 
tentials are given. Since our treatment is for general values of N, the results can be studied 
in the large N limit. A somewhat surprising result is that a free particle incident on an 
A-prong vertex undergoes continuously increased backscattering as the number of prongs 
is increased. 
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1 Introduction 



Various aspects of the solutions of the Schrodinger equation for both scalar and vector 
potentials on a wide variety of networks have been discussed by several authors [1-8]. Re- 
cent investigations have been motivated in part by the considerable interest in mesoscopic 
systems and the experimental observation of persistent currents Q. Most of the results 
correspond to piecewise constant scalar potentials V = Vq 1^, vector potentials A as- 

sociated with uniform magnetic fields Q, or 5— function potentials at the network vertices 
||2|, ^ 1^. In this paper we discuss the bound states and the scattering properties of a single 
particle moving in an arbitrary scalar A^-prong potential. The treatment is kept simple in 
order to clearly show the generalization from the two-prong case, which is just the familiar 
one-dimensional (particle on a line) problem. The properties of the eigenstates of A/^-prong 
potentials are not a priori evident, since this system is in some sense intermediate between 
one and two dimensions. 



The plan of this paper is as follows. In Sec. 2, we define A^-prong potentials and the 
boundary conditions which the wave function must satisfy. The normalization of wave 
functions and their orthogonality properties are also discussed. Sec. 3 contains a discussion 
of bound states for potentials with A'' identical prongs. The general solution is quite easy to 
obtain and is in fact closely related to the solution of a symmetric one-dimensional potential. 
Sec. 4 contains a discussion of bound states in A'^-prong potentials with non-identical prongs. 
Several analytical and numerical solutions are given to illustrate eigenfunction properties. In 
particular, we show that a generalized version of the usual rule about one extra node for each 
higher eigenstate is probably true but the theorem about non-degeneracy in one-dimensional 
problems is not. Since the concept of supersymmetry has yielded many interesting results 
for one-dimensional quantum mechanics 1 10, 11 1, it is natural to investigate what new results 
come from using supersymmetry for A^-prong potentials. This is discussed in Sec. 5. Sec. 
6 illustrates the use of lowest order perturbation theory to obtain bound state energies. 
Scattering is discussed in Sec. 7. For the special case of no potential in any of the A^ prongs, 
very interesting properties are obtained which are quite different from naive expectations. 
Finally, in Sec. 8, we extend the usual discussion of tunneling for a double well potential 
to the case of A^-well potential systems. 



2 Multi-Prong Potentials 

We are addressing the problem of formulating and solving the Schrodinger equation for 
a single particle constrained to remain in a space made up of A^ lines meeting at a vertex 
point (Fig. 1). The prongs are labeled by the indices i (i = 1, 2, • • • , A^) and the position on 
any prong is given by the positive coordinate Xi, with Xi = being the vertex point. The 
potential is fully specified by giving the scalar potentials Vi{xi) on each prong i. The overall 



2 



wave function ip{x) is composed of the individual wave functions ipi{xi) on each prong i: 



iP{x) = {iPi{xi),---,iI)n{xn)]- (1) 

We are using the arrow notation to indicate the appropriate A'^-tuple. The physical require- 
ment of single- valuedness of the wave function at the vertex implies for the component wave 
functions 

V'i(0)=V'2(0) = --- = V'^(0). (2) 

In this paper, we will not be considering J- function potentials at the vertex. Consequently, 
the second vertex condition requires the sum of all derivatives to add up to zero [ffl, 0] : 



N J / 



= 0. (3) 

Xi=0 



This can be readily derived as follows. The time independent Schrodinger equation, reads 
V'^ip = {V — E)ip, taking units in which h = 2m = 1. We assume that the potential 
is finite in the neighborhood of the vertex, and integrate over a small sphere of radius 
e with the vertex as the center. Using the divergence theorem and letting e — > yields 
lime^o Ja ■ ~ 0) which reduces to Eq. for an A^-prong potential. Note that for 
the special case N = 2, one has the standard one-dimensional situation and, taking into 
account the relative directions of xi and X2, Eq. @ is just the statement that the derivative 
dip/dx is continuous across the vertex. 

For the A^-prong system, we call two functions %l){x) and </>(x) orthogonal if 

N 

/ ip*{xi)(l)i{xi)dxi = 0, (4) 

i=i 

where the sum indicates that all prongs are included and each integral runs over the appro- 
priate prong length. The normalization condition, which also includes all prongs, reads 

N 

"Y i;*{xi)'il^i{xi)dxi = l. (5) 
i=i •' 



The eigenfunctions of the Schrodinger equation in a multi-prong potential not only have 
to satisfy the vertex conditions given in Eqs. (§) and (|3|), but also appropriate boundary 
conditions at each prong end. These conditions depend on the energy E and the specific 
behavior of Vi{xi) for the maximum allowed value of Xj. Eigenfunctions corresponding 
to bound state solutions have to vanish at the prong ends. For potentials which reach 
a constant value fast enough at the prong ends, one has standard plane wave solutions 
^^ikxi ^ ^-ikxi^ for the scattering states. 
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3 Bound States: Identical Prongs 



We begin by discussing the interesting special case of N identical prongs, i.e. Vi{x,i) = 
V{xi). Note that this does not imply that the wave functions are identical on all the prongs. 
In fact, Eq. (|3|) implies the wave functions are often different, since their derivatives at the 
vertex have to add up to zero. This statement is the generalization of the two-prong case, 
where one has symmetric and antisymmetric solutions. 

The general bound state solution for N identical prongs is easily obtained from the 
two-prong situation. For two identical prongs, one has the situation schematically shown 
in Fig. 2. We can define the single variable x such that x = xi for positive x, and 
X = —X2 for negative x. Effectively, one is mapping the two prongs onto the real axis, 
— oo < x < oo. Clearly, this is the familiar situation of a symmetric, one-dimensional 
potential. Its eigenstates correspond to even and odd solutions ^^'^\x) at energies En, 
labeled by a quantum number n (n = 0, 1, 2, • • •). 

Theorem: The eigenstates of a potential with N identical prongs can be constructed from 
the eigenfunctions ip^'^\x) of the corresponding symmetric two-prong system, and have 
exactly the same eigenenergies. Explicitly, 

V^(")(a,x) ^ {aiV'(")(xi),a2V'(")(x2),---,a^V^")(x,v)}, (x. > 0). (6) 

where a is a compact notation for the A^-tuple (ai, 02, • • • , cat). For even numbered states 
the boundary conditions at the vertex imply ai = 02 = • • • = a^', for odd numbered states 
one has the constraint J2i — 0) ^'^d there is an {N — 1) fold degeneracy. 

Proof: By construction, the wave function 'ip^"'\a, x) for an eigenenergy En satisfies all the 
boundary conditions at the prong ends. It only remains to show that the vertex conditions, 
Eqs. @ and (|3|) are also satisfied by Eq. (|6|). a) For even eigenfunctions of a two-prong 
potential one has ip'{0) = 0, and in general ^(0) 7^ 0. Eq. @ is therefore trivially satisfied, 
and Eq. (^) implies oi = 02 = • • • = qat. The eigenfunction is therefore determined 
upto one overall normalization constant, and the energy level is non-degenerate, b) For 
odd eigenfunctions of a two-prong potential, one has fp'^O) 7^ 0, and tp{0) = 0. Now, 
Eq. (|2|) is trivially satisfied (since each wave function vanishes at the vertex) and the 
derivative condition Eq. (|3|) implies J2i — 0; claimed. This constraint is not sufficient 
to fully determine the constants Oj, and the levels have an (A^ — 1) fold degeneracy. The 
completeness of states in Eq. ^ seems intuitively reasonable, and can be established from 
the completeness of the eigenfunctions of a symmetric one dimensional potential. 

As a first illustrative example, consider a particle of mass m moving in a harmonic 
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potential with three identical prongs 

= ]^mu?xi. (7) 

Each prong corresponds to the positive half of a harmonic oscillator potential. The ground 
state of this three-prong potential will have the same energy £"0 = ^^a; as the two-prong 
case. The un-normalized eigenfunction is given by 

V'(°)(a,x) = {aV'(0)(xi),aV'(°)(x2),a?/;(0)(x3)}, (8) 

where, ij)'^'^'^ oc e"'^^"/^ is the Gaussian harmonic oscillator ground state wave function. We 
assume il>^^\x) to be normalized to unity over one prong, so that the wave function in 
Eq. (P) is properly normalized. The next level at energy Ei = ^fiuj is doubly degenerate. 
Eigenfunctions at this level are given by 

V'(') (a, x) = {ai^/.(i) {xi),a2^^^^ (xs), as^^'^ (^3)} , (9) 

where, V^"^^ oc xe""^^^/^ is the wave function for the first excited state of the two-prong 
case. There is an infinite number of wave functions ip^^\a,x) that correspond to the same 
energy ^hui. Each of these wave functions is characterized by a set of coefficients obeying 
the constraint J2i — 0- Due to this constraint, we see that all allowed span a two- 
dimensional space, thus leading to a two fold degeneracy at this level. The following two 
functions then provide one choice for an orthonormal basis in this two-dimensional, two-fold 
degenerate space: 

(10) 

This pattern of nondegeneracy/degeneracy keeps repeating as we consider higher energy 
states. All even numbered states are nondegenerate and all odd numbered states have two 
fold degeneracy. The three lowest eigenstates are shown in Fig. 3. The normalization has 
been chosen such that il)'^'^'^* {x)^^^\x)dx = 1. 

Similarly, the eigenfunctions for the potential V{xi) =0, Xi < 1 are shown by the 
dashed lines in Fig. 4. The wave functions are sinusoidal, and the energies are En = 
n^7r^/4 (n = 1, 2, • • •), since they correspond to a symmetric two-prong potential which is 
an infinite square well of width 2. 

From the above examples, it becomes clear that a potential with N identical prongs 
will have alternate nondegenerate and {N — 1) fold degenerate energy levels. Clearly, the 
familiar non-degeneracy property of one-dimensional potentials is maintained for N = 2, 
but not for situations with more prongs. 
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4 Bound States: Non-Identical Prongs 



Let us now consider the general case of an A^-prong system with potential Vi{xi) on 
the i-th prong. Let ij^i^Xi) be the solution of the Schrodinger equation with energy E along 
that prong: 



+ {Vi - E) 



0, 



{i = l,2,---,N). 



Of the two linearly independent solutions to this equation, let (j)i{xi,E) be the one that 
vanishes at the end point of the i-th. prong. This implies ^i{xi,E) = ai<j)i{xi, E), where aj 
is a constant. These Tpi{xi, E), glued together properly so as to satisfy vertex requirements, 
will produce the wave function for the entire domain. From vertex conditions (|2|) and (|^, 
we get 

MO,E)=aiUO,E) = U, (11) 



and 



N 



N 



Y^^[{Q,E)=Y^a,c^[{Q,E) = f). (12) 

i=l i=l 

The energy eigenvalues E are then determined by eliminating the unknowns Oj and U from 
the above set of {N + 1) linear equations. The eigenvalue condition is: 



det 



h{Q,E) 






^2(0,^) 





0;(o,i?) ,^'2(0, i^) 








<I)n{Q,E) 

4>'n{^,e) 



1 
1 

1 





0. 



(13) 



For any given eigenvalue En, Eq. (^) then determines all the constants in terms of one 
unknown constant U . The corresponding bound state wave function is: 

4>l{xi,En) (j)NixN,En) 



^{x,En)=U 



(14) 



0l(O,^„) ' ' (t>N{^,En) 

If desired, the constant U can also be fixed by requiring overall normalization of the wave 
function. 



To illustrate the procedure for determining eigenvalues, consider the example of a par- 
ticle that is free to move inside a domain of three prongs of lengths li, I2 and l^. Since the 
potential along the prongs is zero, the wave functions along them are given by sinusoidal 
functions that vanish at the end points. An eigenfunction of energy E = k"^ has the form: 

'4'{a, x) = {ai sin k{li — xi), 02 sin k{l2 — 2:2), 03 sin k{l^ — X3)} . (15) 

In this case, Eq. (|Tl| ) is: 

ai sin kli = 02 sin kl2 = 03 sin kl^ = U, (16) 
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where U is the common value of the wave function at the origin. The derivative condition 
dH) gives: 

ai cos kli + 02 cos kl2 + 03 cos kl^ = 0. (17) 
The eigenvalue condition ([l3|), when simplified reads 

cos kli sin kl2 sin /cZs + sin kli cos sin kl^ + sin /c/i sin kl2 cos /c/a = 0. (18) 

As a specific case, consider the situation I2 = 0.8 and /i = ^3 = 1. For this case, 
the solutions to Eq. (|l8D are found to be /c = 1.68, vr, 3.61, 5.08, 27r, 7.17, 8.54, Svr, .... The 
solutions k = mir (m = 1,2,...) are a consequence of maintaining partial symmetry by 
taking two prongs to be identical {h = I3 = 1). The corresponding wave functions are 
shown by the solid lines in Fig. 4. [For comparison, we have also plotted the dashed lines 
corresponding to the wave functions for the case of all three identical prongs of length 1]. 
From the figure, it is apparent that the ground state wave function at Eq = 1.68^ has no 
nodes. The first excited state at Ei = vr^ has one node at the vertex point. Note that 
for this state, the wave function in prong 2 is zero. The second excited state is at energy 
E2 = 3.61^ and has two nodes, one in prong 1 and the other in prong 3. This result is 
very suggestive - one expects one extra node to appear for each higher eigenstate, similar 
to the familiar one-dimensional situation. 

We have also studied the variation of the eigenvalues k systematically as the prong length 
I2 is varied. The results are shown in Fig. 5. As discussed above, the solutions k = rmr are 
present, and when I2 is an integral multiple of /i(= ^3), one has the interesting occurrence 
of degeneracy and level crossing. The curves in Fig. 5 are labeled by the number of nodes 
in the wave functions. At any fixed value of I2, the number of nodes increases with energy. 
Note that as I2 0, the eigenvalues become doubly degenerate at A; = mvr, (m = 1, 2, ...). 
This is intuitively clear since the limit /2 ^ forces the wave function to vanish at the 
vertex, thereby effectively breaking the problem into two infinite square well potentials of 
widths h = 1 and ^3 = 1. Finally, in Fig. 6, we take li ^ I3, and plot eigenvalues k. This is 
a completely asymmetric situation, and there is now no degeneracy of energy levels. Again 
note that the number of nodes in the overall wave function increases by one with increasing 
energy levels. 



As a second example, we determine the eigenstates of a three-prong potential Vi{xi) 
juifxf, composed of three harmonic oscillators of different angular frequencies uJi {i 



1, 2, 3), and h = 2m = 1. The wave functions which vanish at Xj ^ 00 are [13 



E 1 

i/jiixi) = aiD,,-{^iXi) , Vi = , (19) 

where Dy is a parabolic cylinder function. The eigenvalue condition [Eq. ([l^)] gives 

V^TI?!,^ (0)1)1,2 (0)-Di.3 (0) -h cyclic permutations = 0. (20) 
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The parabolic cylinder function and its derivative at the origin have simple expressions in 



terms of gamma functions |13]: 



The solutions of the eigenvalue equation are then easily found. For the choice loi = 1.0, uj2 = 
2.0, 1^3 = 3.0, one gets the four lowest eigenstates at energies E = 0.83,1,94,3.29,3.85. 
There is no degeneracy since all w's are different. 

In the above two examples, exact analytic forms for the solutions ipi{xi) were available. 
Even if this is not the case, it is easy to use numerical Runge-Kutta techniques applied to 
each prong. 

5 Supersymmetric Quantum Mechanics on Multi-Pronged 
Domains 

Given any one-dimensional potential V{x), the powerful techniques of supersymmet- 
ric quantum mechanics can be used to generate a partner potential V{x) with the same 
eigenvalues |10, 11]. This property has been extensively used to get a deeper understanding 
of exactly solvable potentials and for generating improved approximation methods (SWKB 
method |14], large method |15|, etc.) for determining eigenvalues. Supersymmetric quan- 
tum mechanics provides an elegant formalism which includes and goes substantially beyond 
the method of factorization previously applied to some potentials fl^. In this section, 
we show how the ideas of supersymmetric quantum mechanics can be applied to a given 
A^-prong system V{x) in order to generate solutions for a new A/'-prong potential V{x). 

As in the previous section, consider a scalar A'^-prong potential V{x) composed of po- 
tentials Vi{xi) on prong i (i = 1, 2, A^). Its eigenvalues and eigenfunctions are then given 
by Eqs. ( p!3| ) and (14) respectively. The un-normalized ground state wave function is 



^ ^"'''°^"i<)(o,i^o)' ' 0S^(o,i^o) r ^ ^ 

This can be used to define the superpotential W{x) whose value on prong i is 



W,{x,) = -^^^. (23) 

It is easy to check that 

Viixi) = Wlixi) - W[{xi) + Eo. (24) 
The supersymmetric partner potential is given by 

VAxi) = W,^{xi) + Wlixi) + Eo. (25) 



8 



The potentials Vi{xi) taken together make up the full iV-prong potential V{x). From su- 
persymmetric quantum mechanics, we know that the solution of the Schrodinger equation 
for potential Vi{xi) and energy E is given by |11| 



(26) 



where Oj are constants. Since (pi{xi,E) vanishes at the prong ends, so do all the (l)i{xi,E). 
At the vertex, one wants 



U0,E)=4>j{0,E) and ^<^^(0,^)=0. 

i 

This gives an eigenvalue condition similar to Eq. ([l^). 



(27) 



det 










^2(0,^) 



UO,E) 








J)n{0,E) 1 
i'j^{0,E) 



0. 



(28) 



In general, the eigenvalues obtained from Eq. ( pS]) will be different from those coming 
from Eq. (|l3|). However, there are two important situations where the same eigenvalues 
result. This happens for the special case of two prongs {N = 2) and for the case of all 
identical prongs. These results can be understood from a physical viewpoint, since N = 2 
is the standard one-dimensional situation treated in supersymmetric quantum mechanics, 
and furthermore, as we have seen in Sec. 3, these are the same eigenvalues for the identical 
prong case. Thus we see the machinery of supersymmetric quantum mechanics can be 
immediately applied to get eigenvalues and eigenfunctions of the partner potential V{x) for 



the identical prong case. The equality of eigenvalues from Eqs. (13) and (|28| ) can also be 
established from a mathematical viewpoint using Eqs. (p^) and (27). 



6 Perturbation Theory 

Having solved a multi-prong problem, it is of interest to see the influence of a small 
perturbation of the potential in one of the prongs, say, on the energy levels. Here, for 
simplicity, we will only deal with an unperturbed system with three identical prongs. Con- 
sider an unperturbed system with potential Vi{xi) = Vo{xi) on each prong. Its completely 
symmetric ground state 
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is described in terms of tp^^\x); the ground state eigenfunction of a two-prong system 
with the same potential Vo(xj) in each prong. We have chosen il^^^^x) to be normalized 
to unity over one prong, so that the wave function in Eq. ( p9[ ) is properly normalized. 
The eigenenergy of this state is also the same as the ground state energy of the two-prong 
system. The first excited state for this three-prong system is doubly degenerate. We choose 
the following two orthonormal states as our base states: 

V^«(ai;x) ^ i^{v,«(xi),0,-V(^)(x3)} , (30) 



and 



^(i)(a2;x)^-^{V'«(xi),-2V'(i)(x2),V'(i)(x3)}. (31) 

Both of these states have exactly the same energy as the first excited state of the two-prong 
system. 

Now, let us include a perturbation Vj that is nonvanishing in only one prong. The 
energy of the system is shifted and the degeneracy of the states of Eqs. ( |30[ ) and ( |3l] ) is 
lifted. Interestingly, only one of the eigenvalues of the first excited state changes, while the 
other one remains the same, thus causing the split. We compute here such shifts for the 
ground state as well the first excited state of the system. For numerical concreteness only, 
in the following example, we use the harmonic oscillator potential for Vq and axf as the 
perturbation potential Vj on prong 1. Thus the unperturbed potential is 

Voixi) = )u?x\ . (32) 

Following the usual route, the shift of the ground state can be computed. The first order 
shift, (JE'o, is given by 

5i?o = 1 dx^^f\x{)Vi^f{x{) = ^-^^ . (33) 

The next energy level has two fold degeneracy, and hence it is necessary to use the 
formalism of degenerate perturbation theory. We compute the matrix elements of the 
perturbing potential V/ using -i/'^^n^l; ^) and V^^-^H^^; ^) as basis vectors. For this situation, 
the matrix Vi is explicitly given by 



^7 h 



where, 7 = ., . Diagonalizing the corresponding Hamiltonian, we find the eigenvalues 

V 27ra;^ 

to be and I7. Hence, at this level degeneracy is lifted and the difference in energy is ^7. 
The two eigenfunctions are given by 

^^«(al;x) - ^^^^){a^-x) = -L{0, (^2), (xg)} (34) 
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and 

i^«(al;f) + ^^{i)(a-2;x) = i= {2V'«(xi), -^^^(xs), V^^^^s)} • (35) 

It is interesting to note that the first state has the same energy as the unperturbed system. 
Thus we see that, since the perturbation was added to just one prong, it does not break the 
symmetry completely and one of the eigenenergies remains unchanged. 



7 Scattering in a Multi-Prong System 

The scattering of an incident plane wave off the vertex in a multi-prong system offers 
scenarios that are different from the usual scattering in one-dimensional problems. As a 
first example, consider a plane wave with energy E moving along Prong 1, incident upon 
the vertex of an A^-prong system with a constant potential y = on all prongs. This is 
a trivial example in the two-prong case leading to 100% transmission, and zero reflection. 
For more than two prongs, it is not a priori intuitively clear whether full transmission will 
occur or not. The wave function on prong 1 is given by 

xpi{xi) = exp(— i/cxi) + r exp(z/cxi), (36) 

where k is the momentum of the incident plane wave {E = k'^). The wave function on all 
the other prongs {i ^ 1) only consists of outgoing waves and is given by 

itiixi) = ti exp{ikxi). (37) 

Imposition of the boundary conditions (|2|) and on these wave functions at the vertex re- 
lates the reflected amplitude r to the transmission amplitudes ti . Their explicit relationship 
is given by 

N 

{l+r) = ti {i = 2,---,N), r-l + J2ti = 0. (38) 

i=2 

The solution is 

2- N 2 

U = - {i = 2r--,N). (39) 



N ' N 

The reflection and transmission coefficients are then given by 

Clearly, one has TZ + J2iL2'^i — 1) the probability current is conserved. For N = 2, 
the reflection coefficient vanishes, as expected. However, for N > 2 there is always a finite 
amount of reflection. For example, for a three-prong potential, the reflection coefficient is 
71 = ^. In fact TZ increases continuously with A^, and approaches unity for large number of 
prongs, A^. (This increase of TZ with A^ is shown in Fig. 7 by the curve labeled 1.) Thus 
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one has the curious result that an incident wave, when given a large choice of scattering 
paths, in fact prefers to be reflected back on its initial prong! 



Now let us consider a more general case where the incident wave is on a prong with 
potential zero, whereas the remaining (A^ — 1) prongs are at constant potential Vq. Let the 
energy of the wave be ii^ > Vq. Let us define k"^ = E, k''^ = E — Vq and the parameter 
C = = ^Jl — The wave function on prong 1 is again given by ^"1(2^1) = exp(— iArxi) + 
r exp{ikxi), and wave functions along remaining (N — 1) prongs are given by ipi{xi) = 
tiexp(ik'xi) (i = 2, • • • ,N). Boundary conditions (§) and @ applied at the vertex yield 

(l + r) =t2 = t3 = ■■■ = tN = t , - 1 + r + (TV - 1)C = 0. (40) 

Solving these equations, we find r = and t = jqi^^^^r^j. One can show that r 

and t obey 1 = r-^ + (^(iV — l)t^, as required by conservation of probability. The reflection 
coefficient TZ, which is now given by 

2 



7^ 



i-C(iV-i) 



i + C(iV- 1) 



has a rather interesting behavior as a function of ^ and A^. In Fig. 7, we have plotted TZ 
as a function of N for several values of TZ vanishes for C = jj^^. Thus for some special 
values of we have complete transmission through the vertex. As N increases, TZ slowly 
approaches unity for all 

As yet another interesting example we consider scattering from a T-stub with open 
ended prongs and one prong of finite length I (the stub). We consider a case where potentials 
along all the prongs are zero. A plane wave of momentum k (energy E = k"^) incident upon 
the vertex from prong 1, sets up stationary waves in the stub, and outgoing plane waves 
in the remaining — 1 prongs. The wave function on prong 1 is given by exp(— i/cxi) + 
rexp(iA;xi). Along the stub, the wave function is given by Asva.k{l — Xg), where Xs is the 
coordinate along the stub and A is a constant. Wave functions along all other prongs are 
given by texp{ikxj), {j = 2, ■ ■ ■ , N). From boundary conditions (§) and (^), we get: 

l + r = Asinkl = t ; ik{r - 1) + ik{N - l)t - Ak cos kl = 0. 

Solving them, we get 

-{N -2) - cot kl 2 



N + icotkl ' N + icotkl' 

|2 _ cot '^kl+{N -2)'^ 



(41) 



The reflection coefficient is given by 7^ = |r| = — cot '^ki+N'^ • Clearly, from the above 

expression, one gets |rp + (A^ — l)|tp = 1. Interestingly, we find that for kl = rnr, i.e. 
2 2 

energy E = "^J , the reflection coefficient 7^ — > 1. At these energies, a standing wave is set 
up jointly in the stub and in prong 1, and the particle never enters the other prongs. 
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Let us generalize the discussion of scattering to an A'^-prong vertex with identical po- 
tentials V{xi) on all prongs. Define two linearly independent solutions /(xj) and g{xi) of 
the Schrodinger equation using asymptotic boundary conditions as x, ^ oo: 

f{xi) exp{-ikxi), g{xi) exp{ikxi). (42) 

The vertex conditions are /(O) + r g{0) = t g{0), and /'(O) + r g'{0) + {N -l)t g'{0) = 0. 
The solution is 

_ i/'(o) (iv-i)/(o) ^ if'{o)^i_m^ ^^^^ 



Ng'iO) N 5(0)' Ng'iO) N g{0) 

Conservation of probability requires |rp + (A^ — l)|tp = 1 for all values of N. This 



is satisfied provided = = 1. Here again, we see that as the number of prongs 

N —>■ oo, the reflection coefficient TZ = |rp approaches unity. 

Having analyzed scattering on a general identical A^-prong domain, we now investigate 
whether we can extract any further information using supersymmetry. In one dimensional 
quantum mechanics, supersymmetry relates reflection and transmission coefficients of one 
potential with those of its supersymmetric partner potential! 12 1. We find that this relation- 
ship holds also for the case of N identical prongs. 

We have just seen that for the potential V{xi) on all prongs, the scattering is described 
by Eq. (43). The partner potential is given by Eq. (^). Let / and g be the solutions of the 
partner potential which behave like /(xj) — > exp(— ifcxj), g{xi) exp(ik 
Functions / and g are given by 

where Woo is the value of the superpotential at infinity. From these expressions, we can 
show that = (^) f|) , and = (^) The reflection and transmission 
amplitudes for the potential V{xi) are given by 

fWoo-ik\f 1/(0) {N-l)f'{0)\ ~_ / Woo-ik \ / 1 /(O) ^ 1 f(0) 



\Woo+ikJ\ N g{0) N g'{0)J \Woo + ik J \ N g{0) Ng'{0), 

Note, for the special case of two-prongs, these equations agree with reference jl^. We 
flnd the following relationship among r, t and their partners r and t: 

Woe -ik\^ r Woo-ik -^ 



Moo + ik) ' \Woo + ik, 

Note that |fp = |rp, which implies that for the identical prong case, the reflection coeffi- 
cients are the same for any potential and its supersymmetric partner. 

Now, using an explicit example, we provide a concrete demonstration of the usefulness 
of the machinery that we just developed. For the sake of simplicity, we shall work with 
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the example of a free particle in an A'^-prong system for which reflection and transmission 
coefficients are given by Eq. (]39|). The superpotential for this problem is given by W{x) = 
tanhx. This superpotential generates two distinct potentials related by super symmetry. 
They are V{x) = W'^{x) - W'{x) = 1 -2sedi'^x and V{x) = W'^{x) + W'{x) = 1. The first 
potential holds one bound state at energy E = 0. 

Amplitudes r and t for the free particle system are and ^, respectively. Amplitudes 
r and t for the potential V{x) = 1 — 2sech then given by 



1+ik 
1-ik 



1 /(O) N - 1 /'(O) 



Ng{0) N g'{0) 



l + ik\ /N -2\ /l + ik\ ( 2 



\-ik \ N ' \\-ik \N 



Thus, knowing the reflection and transmission amplitudes f and i for the free particle 
system, supersymmetry allows us to determine amplitudes r and t for the rather non-trivial 
potential V{x) = 1 — 2sech^x. Note that coefficients TZ = |fp and T = |ip are equal to 
coefficients for the partner potential, TZ = |rp and T = |tp respectively. 



8 Tunneling in a Multi-Prong System 

Tunneling is another sector of great interest that shows a marked difference from the 
two-prong case. Consider a system with N identical prongs, each one having one minimum 
(well-like structure, similar to the two-prong case shown in Fig. 2). In this section, we look 
at the tunneling of a localized wave packet from one prong to the others. For simplicity, 
we start with a three-prong system. The generalization to higher number of prongs is then 
straightforward. Let denote a wave packet localized in the i-th prong. Such a packet 

can be approximated by the following linear combination of eigenstates of the system: 

= i + V(^)(a-;,x)) , (45) 

where 1 = (1, 1, 1) and denotes a three-tuple that has a 2 in the i-th position and —1 in 
the rest. For example. Si is given by the three-tuple (2, —1, —1). 

The amplitude for the localized state $j(x) to tunnel out of the prong 1 is given by: 

A{t) = {iPf\exp{-iHt)m (46) 

where 

V'i = ^>i(x) and V/ = ^^^= (a2^2(^) +a3^*3(^)) ; 



«2 + 03 = 1- Substituting explicit expressions for ^l;f and ^/^j in Eq. ([46[), we get 
^„,„3 = , ^ / dx [a2$2(x) + a3$3(x)] fe-^^" V^°Hl,^') + e'^^^ V^'^(«"l, i') 
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9^/ai + a§ 



dx 



a2tp'^^^ {a2,x) + a3^(^) (al, x) + i/j^^^ (1, x) 



a2di ■ (12 + "sOi • + 3e"^^* , 5E = Ei - Eq 



where we have used / ■(/''•"^H'^! ('^'2; s;) = 01-02 = Z]i(o'i)i(a'2)i. Substituting ol-o's 

oi • a'2 = —3, we get 



A, 



Oi2 Oi3 



The probabihty of tunnehng is then given by 



-1 + e 



iSEt 



a2 0:3 



= \A, 



■Oi2 CX3 I 



/6Et 



sm 



The maximum value for this tunnehng probabihty is |, and that occurs when a2 = = ^■ 
It is worth noting here that unhke the two-prong case, the wave packet never completely 
goes out of the prong 1. 

When we generalize this situation to A^-prong case, we find that the probability ampli- 
tude for a wave packet of tunneling out of prong 1 is given by 



Aa 

where a = {02, as, • • • , (Xn} with 



2ie-'^E"+^^/^hm{SEt/2) 

giving the amplitude for finding the packet in 



This amplitude is maximum for the symmetric case for which 



the i-th prong 

all ai are equal. In this case the probability of tunneling is given by 

Thus the probability of tunneling decreases when more alternatives are available, and goes 
to zero as the number of prongs becomes very large. This is similar to the scattering 
situation discussed in Sec. 7. 
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Figure Captions 



Figure 1: Schematic diagram of an A^-prong potential. The position on any prong i is 
specified by the positive coordinate rcj, with Xj = being the vertex. The angles between 
prongs play no role in any computations. 

Figure 2: Example of a potential with two identical prongs. The positive variables xi,X2 
together span the real x-axis (— oo < x < oo). 

Figure 3: The three lowest eigenfunctions [Eqs. (|8|), @ and ([lO|)] for a three-prong 
harmonic oscillator potential V{xi) = ^muj'^xf. The ground state if)^^^ is non-degenerate, 
whereas the first excited state 'ip^^^ is doubly degenerate. 

Figure 4: The three lowest energy eigenstates (un-normalized) for a three-prong free par- 
ticle potential. The solid lines correspond to h = h = 1, h = 0.8 . The dashed lines 
correspond to the identical prong case li = I2 = I3 = 1 . 

Figure 5: The eigenvalues k (energy E = k"^) for a three-prong free particle potential as 
a function of I2, the length of prong 2. The lengths of the other two prongs are kept fixed 
and equal li = = 1 . Note the pattern of degeneracy and level crossings. The curves are 
labeled by the number of nodes in the wave function (in all prongs). 

Figure 6: The eigenvalues k (energy E = k'^) for a three-prong free particle potential as 
a function of I2, the length of prong 2. The lengths of the other two prongs are kept fixed 
and unequal li = 1 , l^ = \pl. There is now no degeneracy. The curves are labeled by the 
number of nodes in the wave function (in all prongs). 

Figure 7: A plot of the reflection coefficient 7^ versus the number of prongs N . The 
incident wave with energy E is on prong 1 with zero potential. All other prongs are taken 
to be at a constant potential Vq. The curves are labeled by the parameter = wl — 
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